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PACS 45 . 70 . -n - Granular systems 

PACS 72 . 80 . -r - Conductivity of specific materials 

PACS 81.05.Rm - Porous materials; granular materials 

Abstract. - A DC electrical current is injected through a chain of metallic beads. The electrical 
resistance of each bead-bead contacts is measured. At low current, the distribution of these 
resistances is large and log- normal. At high enough current, the resistance distribution becomes 
sharp and Gaussian due to the creation of microweldings between some beads. The action of 
nearby electromagnetic waves (sparks) on the electrical conductivity of the chain is also studied. 
The spark effect is to lower the resistance values of the more resistive contacts, the best conductive 
ones remaining unaffected by the spark production. The spark is able to induce through the chain 
a current enough to create microweldings between some beads. This explains why the electrical 
resistance of a granular medium is so sensitive to the electromagnetic waves produced in its vicinity. 



The electrical resistance of a granular assembly is very 
sensitive to a large variety of external perturbation. The 
global electrical resistance can be indeed modified by a 
mechanical shock or stress, by a thermal dilatation, by 
aging [1], by applying an electrical current [2,3], or by 
l!)roducing electromagnetic perturbation in its vicinity [4] . 
The two last sources of perturbation are the most unex- 
I^ected. The relation between the voltage and the current 
injected through a granular material has been debated for 
a long time. It has been shown at the end of the XlXth 
century that a huge decrease of resistance occurs when a 
current is injected through metallic fillings [5]. The re- 
sistance drops over several order of magnitude when the 
current reaches a given threshold. Almost during the same 
period, E. Branly discovered that an electromagnetic wave 
(e.^., spark production in the air) is able to modify the 
electrical resistance of a granular heap at distance [4]. 
This remarkable phenomenon is at the origin of the de- 
velopment of the wireless transmission. These problems 
have been recently revisited because the mechanisms are 
still not completely elucidated [6]. Moreover the electri- 
cal properties could be a smart way to probe the internal 
structure of the mechanical arches through a granular pile. 
When an electrical current is injected through a pile, it 
percolates according to the least resistance pathway which 
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Fig. 1: Schematic view of the experimental setup. The beads (8 
mm of diameter) are placed into a groove and electrodes have 
been soldered on each bead. A current can be injected through 
the chain via the beads located at the extremities. Sparks can 
be also produced between two electrodes separated by 2 mm 
and located at a distance d from the chain. 



has a topological dimension of 1. One of the challenges is 
to determine the exact role of the network compared to 
the role of one single contact with respect to the imposed 
perturbation. 
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Fig. 2: Typical current- voltage characteristics for 4 single contacts within the bead chain (different line styles), (a) Effect of 
a strong applied current : / is first increased from 10 /iA to 1 A, then is decreased to 10 /iA. (b) Effect of sparks: / is first 
increased from 10 /iA to 1 mA, then sparks are produced at 1 m from the chain. The influence of the sparks is shown by the 
discontinuity in two curves. Afterwards, / is increased up to 1 A before being set back to 10 /iA. The single (resp. double) 
arrows denote an increase (resp. decrease) of the current. All the 4 curves collapse on each other when the current is decreased. 
The horizontal dashed line shows the saturation voltage V* :^ 0.4 V (see text). 



The relation between DC current and voltage has been 
described for a one-dimensional (ID) chain of metallic 
beads [2], in a 2D configuration [7,8] or in a 3D pack- 
ing [9, 10]. As reported in these works, the electrical 
properties of the bead assembly strongly depend on the 
electrical history of the granular pile. The voltage is 
not univocally determined by the current because of irre- 
versible processes such as microwelding occuring between 
the beads [2]. Some works have also revisited the influ- 
ence of sparks on the electrical resistance of a 3D packing 
of lead beads (Branly effect) [11] and theoretically [12]. 

The aim of this paper is to study the effect of either 
a high current or an electromagnetic perturbation on the 
resistance of one single bead-bead contact. By performing 
the experiments several times, a large number of contacts 
will be considered in order to establish the distribution 
of the resistances before and after the application of the 
perturbation. The comparison of both situations will allow 
to determine the behavior of one single contact. 

The experimental setup is shown in Fig. 1. Seventeen 
stainless steel beads (8 mm of diameter) are placed in a 
linear groove dug in a nylon block. A screw allows to com- 
press the chain of beads to about 100 N. Electrodes are 
soldered on each bead. Two more wires are soldered on 
the extreme beads in order to inject the current. A sta- 
ble current source (Keithley K2400) is used for this pur- 
pose. A Keithley 2700 multimeter with a multi-channel 
card is used to determine the voltage between the succes- 
sive beads. That ensures a 4- wire measurement for the 
resistances of each contact. Rhumkorf coils are used to 
produce sparks at different distances d from the bead chain 



(d = 0.1 m to 2.2 m). The length of the sparks has been 
fixed to 2 mm, and its duration is roughly 500 ms. This 
ensures that the sparks are produced as soon as the coils 
are switched. In the original Branly 's experiment, an an- 
tenna was fixed to the emitter (Rhumkorff coils) and to 
the granular medium. The antenna allows for the ampli- 
fication of the electromagnetic effects. We decided not to 
use the antenna in order to prevent the masking of any 
effects due to the influence of a low power electromagnetic 
wave. Since the electrical properties of the granular ma- 
terials are very sensitive to their electrical history, after 
each experiment, the system is reset: the pressure on the 
beads is released and the beads are separated from each 
other. Any possible microweldings between beads are then 
broken. The measurements are performed up to 30 times, 
leading to roughly 500 measurements that ensure enough 
statistics. 

The voltages Vi between the bead number i and i + 1 
are measured with respect to the injected current / (see 
Fig. 1). The voltage-current characteristics for 4 different 
bead-bead contacts are displayed in Fig. 2a with differ- 
ent line styles. The initial resistance for each contact 
is defined as the resistance at low current before any ir- 
reversible processes occur (e.^., high current applied or 
electromagnetic wave production). Rq are then extracted 
from each linear fit of the Vi{I) curves between 10 /iA 
< / < 1 mA. The system is reset about 500 times to ob- 
tain enough statistics for the Rq distribution. As shown 
in Fig. 3a, this distribution is found to be very broad over 
4 decades, and is well fitted by a log- normal distribution 
(see •-symbols). A log- normal resistance distribution re- 
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Fig. 3: (a) Probability density functions (PDF) of the resistance Ro at low current (•), the resistance Rs 
a current of 1 A (a), and the resistance Rsp after the production of sparks at 1 m from the chain (■). 
the log- normal distribution fit of i?o (— ) and Rsp 



while a Gaussian is fitted to Rs 



after the injection of 
The curves represent 
The distribution of this latter 



has been plotted on a linear scale in the inset, (b) The cumulative distribution functions (CDF) of Ro 



Rsc ( ), and Rs 



(• • •) are shown on the same semi- log plot. The impact of sparks on the distribution is clearly visible between the Ro (- 
Rsp (• • •) CDF's: only the resistance values of the more resistive contacts have been decreased. 
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fleets the inhomogeneity of the oxide layer on the surface 
of each bead [8]. The cumulative distribution function of 
Ro is plotted in the Fig. 3b (solid line). The mean value 
(Ro) = 38 1^, and its standard deviation is about 155 1^. 
This log-normal distribution will be compared afterwards 
to the one obtained when a constraint is imposed on the 
chain. Two types of constraints can be imposed. Either 
a large current ~ 1 A is injected through the beads, or 
electrical sparks are produced in the vicinity of the chain. 

First, we focus on the effect of a large applied current 
on the voltage-current characteristics. As shown in Fig. 
2a, the current is first increased from 10 /iA to 1 A (sin- 
gle arrows), and then decreased (double arrows). A fine 
measurement of the voltage-current characteristic allows 
to extract several regimes. As shown in Ref. [2,8], when 
the current is increased, three different regimes occur: a 
linear one, followed by a nonlinear part, and then a sat- 
uration regime (see Fig. 2a). When this latter regime 
is reached, the voltage between two successive beads can 
not exceed the saturation voltage V* :^ 0.4 V, and thus 
remains constant when / is further increased [2]. This is 
due to an electro-thermal regulation within the contact. 
In steady state conditions, the temperature of the contact 
can be expressed as [2] 

Tl-TS = ^ (1) 

where V is the voltage across a contact, is the maxi- 
mum temperature reached at the bead-bead contact, and 
To = 290K is the temperature far from the contact re- 
gion, L = 2.45 X 10~^ V^/K^ being the Lorentz constant. 
The contact geometry and the material characteristics do 



not appear in Eq.(l) because both the electrical resistiv- 
ity, p(T), and the thermal conductivity, A(T), are due to 
the conduction electrons, which leads to a linear temper- 
ature dependence Xp = LT, known as the Wiedemann- 
Franz law [13] (see Ref. [2,14] for details). Moreover, the 
size of the micro-contact being much lower than the bead 
size and much larger than the electron mean free path (10 
nm), Eq.(l) holds in a large range of contact size. From 
Eq.(l), a low voltage near 0.4 V increases the temperature 
to about 1370 K at the center of the contact between two 
beads [2]. At such a high temperature, the micro-contacts 
between beads melt leading to microweldings. Since the 
temperature cannot exceed the melting value, it forces the 
voltage, from Eq. (1), to be a constant even when / is 
further increased, leading to a decrease in the resistance 
Rsc = V* / 1 (see the plateau in Fig. 2a). Rsc denotes the 
value of the resistance of a bead-bead contact once the 
saturation regime is reached (that is as soon as a microw- 
elding occurs). When the current / is decreased, the resis- 
tance of the contact then remains equal to Rsc (see double 
arrows in Fig. 2a). The distribution of the contact resis- 
tances Rsc after the passage of a current of 1 A is shown 
in Fig. 3a (see A-symbols). The distribution is found to 
be very narrow, and is roughly fitted by a Gaussian (see 
dashed line in Fig. 3a and in the inset). The mean value is 
(Rsc) =0.29 Since the distribution is very narrow, all 
the contacts can be viewed as electrically equivalent. The 
large current has generated microweldings between some 
beads and has thus erased their initially high resistive val- 
ues. A finer description of the distribution is shown in the 
Fig. 3b by looking at the cumulative distribution function 
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Fig. 4: Logarithm of the ratio between the resistances after, 
i?sp, and before, i?o, the spark production (at various distances 
d from the chain) as a function oi Rq. d — ^.1 (•), 0.25 (0), 1 
(a), 1.4 (v) and 2.2 (□) m. Linear fits of the non-zero values 
of \og{RQ / Rsp) is displayed (— ). Inset: as a function of d. 
Power law fit oc c?^"^ is shown (— ) as a guide for the eyes. 

(CDF) of Rsc The distribution appears slightly asymmet- 
rical due to a limitation at low resistance (the conductivity 
of the steel being finite). 

The intensity of the current has a different effect accord- 
ing to the initial resistance value of a single contact. A 
current is qualified as large once a single contact voltage 
reaches the saturation regime (see Fig. 2a). Once such 
a large current is injected through the chain, the bead- 
bead contact resistances can be split into three groups: 
(i) the highly resistive contacts (with initial resistances 
Rq » V*//*) become much more conducting (their re- 
sistances drop to Rsc ^ V*//*), (ii) those with initial 
resistances Rq « V* /I* are not modified by the current, 
and (iii) the contacts with intermediate initial resistance 
are more or less affected by the current. 

Let us now focus on the influence of sparks on the elec- 
trical properties of the chain. Figure 2b shows typical 
voltage-current characteristics of 4 different bead-bead 
contacts with very different initial resistances Rq (see dif- 
ferent style lines). A current is first applied through the 
chain, and is increased from 10 /iA to 1 mA (single ar- 
rows). Then, sparks are produced using the Rhumkorff 
coils at a given distance d from the chain of beads. This 
leads to a voltage drop of some contacts, the correspond- 
ing Vi - I curves being thus discontinuous. Then, the 
current is further increased from 1 mA to 1 A. The con- 
tact resistances measured just after the spark production 
are named Rsp. Finally, the current is decreased back to 
10 /iA (double arrow). The system is reset several times 
to obtain enough statistics for Rsp . Two kinds of behav- 



iors appear. For inital high resistive contacts (large Rq)^ 
the sparks provoke a sharp voltage drop that betrays a 
drop of resistance (see the solid and dashed lines in Fig. 
2b). The lower contact resistances remain unaffected by 
the spark production (see the dotted lines in Fig. 2b). 
The contact resistances that are influenced by the sparks 
drop to roughly the same value denoted R^^ :^ 2 - 4 1^ (see 
Fig. 2b) and remain stable up to saturation. Note that 
the only resistances much larger than R^^ are affected by 
the sparks. The distribution of the resistance after sparks, 
Rsp^ is displayed in Fig. 3a (see B-symbols). This distri- 
bution is narrower than the initial one, but remains larger 
than the distribution of the resistance after the passage of 
a 1 A current. The Fig. 3b shows the cumulative distri- 
bution function of Rsp well fitted by a log-normal (dotted 
line). Note that the CDF for Rsp for values less than R^^ 
remains the same as the Rq one (see the identical part of 
the solid and dotted lines in Fig. 3b). The arithmetic mean 
of Rsp is equal to 3.8 1^, about 10 times less than {Rq). 

Let us now sum up all of our results. The injection 
of a large current through the chain has the same qual- 
itative effect on the chain conductivity as that by the 
production of sparks in its vicinity. This means that 
sparks can induce a current within the chain enough to 
create microwelding between some beads. One can esti- 
mate the current induced by the electromagnetic waves as 
hnd ^ y /Rip ~ 0.1 A. From Fig. 2a, such current of 
0.1 A is indeed enough to generate microwelding between 
some beads. 

The role of the distance d between the spark emitter 
and the beads on the chain conductivity is now examined. 
Several experiments are performed for different distances. 
Figure 4 shows the logarithm of the ratio between the re- 
sistance after sparks, Rsp^ and the initial resistance, Rq^ 
as a function of the logarithm of Rq. Thus, when the resis- 
tance is not affected by sparks, log(i?o/^sp) equals zero. 
As said above, Rlp{d) is the lowest resistance that sparks 
are able to change. Thus, for a fixed \og{Ro/Rsp) be- 
comes non-zero only for initial resistances Rq > Rlp{d) 
(see Fig. 4). Rlp{d) are then extracted from the linear 
fits of the non-zero data in Fig. 4. The intersections of 
the solid lines with the x-axis give the values of R*p{d). 
The inset of Fig. 4 show the log-log plot of R^p as a 
function of the distance d. i?*^ roughly shows a power 
law dependence on d with an exponent of 1.2. As men- 
tioned above, one can estimate the order of magnitude 
of the current induced by the sparks within the beads as 
hnd{d) ^ V*/Rlp{d) (X 1/d^-'^. This seems to be close to 
the 1/d power-law expected for the electromagnetic waves 
in far field. However, since the bandwidth frequency of 
the emitted waves is unknown, this distance dependence 
deserves further works. When the distance increases from 
(i = 0.1mto2.2m, lind is found to decrease from 0.87 
A to 0.02 A, which remains large enough to produce mi- 
croweldings between beads. 

The infuence of either a high DC current or an elec- 
tromagnetic perturbation on the electrical properties of 



a chain of beads has been studied. The distribution of 
the bead-bead resistances before the perturbation is a log- 
normal over 4 decades. Applying a high current trans- 
forms this distribution to a narrow Gaussian owing to 
the creation of microwelding between the contacts. When 
sparks are produced in the chain vicinity, only bead-bead 
resistances larger than a threshold value i?*^ are affected 
by the electromagnetic waves, i?*^ is roughly proportional 
to the distance between the chain and the spark emitter. 
Spark emission acts as a DC current which intensity in- 
versely depends on the distance between the sparks and 
the contact. This induced current is enough to create mi- 
croweldings between some contacts. Generally, a granular 
packing has a huge contact number. Highly resistive con- 
tacts are thus likely. Since only the largest resistances are 
influenced by sparks, this explains why the conductivity of 
a granular network is so sensitive to the action of nearby 
electromagnetic waves. 
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